A Kinetic Model for Grain Growth 



R. Henseler* M. Herrmann^^ B. Niethammer^ J. J. L. Velazquez^ 

July 22, 2008 



Abstract 

We provide a well-posedness analysis of a kinetic model for grain growth introduced by Fradkov 
which is based on the von Neumann-Mullins law. The model consists of an infinite number of trans- 
port equations with a tri-diagonal coupling modelling topological changes in the grain configuration. 
Self-consistency of this kinetic model is achieved by introducing a coupling weight which leads to a 
nonlinear and nonlocal system of equations. 

We prove existence of solutions by approximation with finite dimensional systems. Key ingredients 
in passing to the limit are suitable super-solutions, a bound from below on the total mass, and a 
tightness estimate which ensures that no mass is transported to infinity in finite time. 

Keywords: grain growth, kinetic model, infinite-dimensional system 
AMS Subject Classification: 35F25, 35R15, 74A50 

1 Introduction 

Many technologically useful materials are poly-crystalline aggregates, composed of a huge number of 
crystallites, called grains, separated by so-called grain boundaries. Typically such materials undergo an 
aging process leading to coarsening of the grain structure and therefore inducing changes in mechanical, 
electrical, optical, and magnetic properties of the material. For details we refer to the review articles by 
Fradkov and Udler |FU94| and by Thompson |Tho01| . 

Different approaches for modelling grain growth in two space dimensions are established in the lit- 
erature. In Monte-Carlo models, compare [ASGS84"al IASGS84b) . the kinetics of the boundary motion 
are simulated by employing a Monte-Carlo technique for moving discrete lattice points describing the 
microstructure. An attractive feature of this model is the simple handling of topological events like grain 
boundary flipping and grain disappearance. 

Boundary tracking models based on partial differential equations as discussed in f KLOU IMNT04) 
offer an alternative to Monte-Carlo models as they deal with quantities of lower dimension. They can 
be further reduced to so-called vertex models where movement of grain boundaries is projected onto 
the triple-junctions, see |KNN89|, IHNO03] . In both cases, however, topological changes require extra 
treatment. 

In the sequel we focus on a kinetic mean-field type models |Fra88( IMar87|, |Fly93| based on the von 
Neumann-Mullins law. Such models consider time-dependent distribution functions for the grain areas 
and the number of sides per grain. Grain areas change according to the von Neumann-Mullins law, 
topological changes are modelled by collision-type operators. Fradkov [Fra88| was the first to develop a 
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model of this type which consists of an infinite-dimensional system of transport equations with a nonlocal 
weight, making the equations nonlinear. In this article we establish a rigorous well-posedness theory for 
this model. 

2 The model 

In this section we present a derivation of Fradkov's kinetic model for grain growth which is based on the 
von Neumann-Mullins law for the change of grain areas and the so-called 'gas ' approximation to describe 
topological changes in a 2D network of grains. Starting point is the isotropic mean-curvature flow for 
the grain boundaries with equilibrium of forces at triple junctions. 

2.1 Networks of grains with triple junctions 

Motion by mean curvature and equilibrium of forces at triple junctions Mean curvature flow 
coupled with equilibrium of forces at triple junctions is a widely accepted model for two-dimensional grain 
growth |BR931 [KL01|, [M!NT04j . For simplicity our objects are 1-periodic spatial networks whose curves 
meet in triple junctions (Fig. [T|). We restrict ourselves in the following to the case of isotropic surface 
energies, such that the curves move according to the isotropic mean curvature flow. Moreover, we assume 
that the mobility of the triple junctions is infinite compared to the mobility of the grain boundaries, and 
this leads to equilibrium of forces at triple junctions. In the isotropic case this condition, also known 
as Herring condition^ just means that the curves meet in an angle of 27r/3. The Herring condition also 
arises as the natural boundary condition in the interpretation of the mean curvature flow of networks as 
L^-gradient flow of the surface energy, see |TC94[ IHROOSj . 




Figure 1: Cartoon of a 2D network of grains with triple junctions indicated by bullets. The Herring 
condition implies all angles to equal 27r/3. 

Von Neumann— Mullins law Under the assumptions stated above (isotropic surface energy, equal 
mobility of grain boundaries, and infinite mobility of triple junctions), one can derive a law of motion for 
the area a (t) of a single grain with n edges [Mul56j , known as the von Neumann-Mullins law: 

±a{t)^Ma^in-6) (1) 

Here M denotes the mobility of the grain boundaries and a the surface tension. The proof uses a direct 
geometric computation involving motion by curvature of the grain boundaries and the prescribed jumps 
of the outer normal by 27r/3 at triple junctions. 

The von Neumann-Mullins law implies that grains with less than six edges shrink, those with more 
than six grow, and such with exactly six edges retain their area (possibly not their shape). 

Topological changes The evolution by mean curvature is well-defined until two vertices on a grain 
boundary collide, after which topological rearrangements may take place. This happens when either 
an edge or a whole grain vanishes. In the first case an unstable fourfold vertex is produced, which 



2 



immediately splits up again, usually in such a way that two new vertices are connected by a new edge. 
In this case, two neighbouring grains decrease their topological class (i.e. the number of edges), whereas 




Figure 2: Neighbour switching 

the two other grains increase it (Fig. The second case causing topological rearrangements is grain 
vanishing. Each grain vanishing is accompanied by disappearance of two vertices and three edges. Due 
to the von Neumann-Mullins law we only take grains with topological class 2 < n < 5 into account. 
Grains with n = 2 and n = 3 vanish in a single possible way. For n = 4 we observe two topologically 




Figure 3: Grain vanishing 



distinguishable possibilities and for n — 5 even five possible local configurations (Fig. For further 
details on the resulting topological classes we refer to the review article by Fradkov and Udler |FU94| . 

It is unclear by which mechanism a specific topological configuration is selected within switching or 
after vanishing events. A natural idea is to compute all possible local configurations and select the one 
that minimises energy locally in the best way, see |HNO03| . 

2.2 Kinetic model for networks of grains 

Our next aim is to derive the kinetic model for large networks of grains with triple junctions. This kinetic 
model comprises the same essential features as the gradient flow dynamics, but differs in some aspects. 

One— particle distribution Following [FraSS) lFU94j we introduce a number density /„ (a, t) that mea- 
sures the number of grains with topological class n > 2 and area a > at time t > 0. Using the von 
Neumann-Mullins law ([T]) we can describe the evolution of / by transport equations 

5J„(a,t) + (n-6)aa/„(a,<) = (2) 

as long as no topological rearrangements take place. Furthermore we choose the following boundary 
conditions 

/„(0,t) = for n>6, 

ensuring that no additional mass is transported from the negative half-axis to the positive one. This 
means no additional grains can be created. 
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'Collision' operator To model topological changes we introduce a collision term {J f)n on the r.h.s. 
of ^ which couples the equations for different topological classes. We define topological fluxes 77+ and 
•q~ denoting the flux from class ri to n + f and from ri to n — 1, respectively, so that 

{J f)n = Vi-l + Vn+l - Vn - Vn ■ 

Next we state the so-called 'gas ' approximation of the collision term by Fradkov |Fra88| that takes into 
account the following transitions between topological classes. 

• Switching events cause simultaneously both a transition from n to {n — 1) (grains that contained 
the vanishing edge) and a transition from n to (n + 1) (grains that contain the created edge). 

• Vanishing of a neighbouring grain corresponds to a transition from n to (n — 1). 

The restriction to these elementary events for topology changes is a simplification of the real dynamics 
in networks of grains. In particular, we ignore that the topological class of a grain is lowered by two if 
the neighbouring annihilated grain was a lens with topological class n = 2. Fradkov and Udler argue 
|FU94j that such an event takes place only very rarely as the number of lenses itself is already very small. 
In addition, we ignore that the annihilation of a grain with topology class n = 5 causes other grains to 
increase their topology class, compare Figured 



Transition rates We make a further simplification by assuming constitutive relations for the transitions 
rates, that in particular ignore all neighbour correlations. More precisely, in what follows we assume that 
the topological fluxes are given by 

?7+=r/3n/„, 7?^ = r(/3+ l)n/„, (3) 

where the coupling weight F describes the intensity of topological changes and depends in a self-consistent 
way on the complete state of the system, see ^ below. The parameter /? in ([3]) measures the ratio 
between switching events and vanishing events. In principle one could allow for arbitrary values of (3 but 
our analysis requires < /3 < 2. Moreover, [FUK88| provide numerical evidence that experimental data 
correspond to 0.45 < < 0.6, that means 0.8 < /3 < 1.5. 

According to Q the collision terms are given by Jf = ^{f)Jf with 

(J/)2=3(/3+l)/3-2/?/2, 

{Jf)^ = iP+l)in+l)fn+i-i2l3 + l)nU + P{n-l)U_i for n > 2. 
Notice that this definition ensures the zero balance property, that is 

(Jf)^^ (a, t) = for all a, t > 0. (5) 

n>2 

This identity reflects that for each ao > the number of grains with area oq does not change due to 
neither switching nor vanishing events. 



Evolution equation The coupling weight F which makes the equations nonlinear (and nonlocal in the 
grain area variable a) is chosen as 

r(/) = rAr(/)/Fz,(/) 

with 

00 00 
riv(/(t)) = E - 6)' fn (0, t), Toifit)) = Y,n f fn (a, t) da - 2 (/3 + 1) / h (a, t) da. (6) 

n>2 n>2 n n 
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We will see in ^'2.'6\ below, that this choice of T guarantees consistency as the solution satisfies the 
polyhedral formula, compare and conserves the total area covered by the grains under the evolution. 



The kinetic model we consider in this paper is thus given by 

dtfn (a, t) + {n- 6) dafn {a, i) = T (/ (t)) (Jf)^ (a, t). (7) 

These equations ([7]) are basically the same as in the work of Fradkov |Fra88[ [FUK88[ [FU94j . The coupling 
term ( J/)2 differs and we do not neglect J f2 da within To (/)• We refer moreover to [BKLT06] . which 
presents some formal analysis and numerical simulations for a similar kinetic model. This model relies 
on different expressions for rj^ but comprises the same essential features as ([7]). 

Finally, we mention the following implication of Definition ([6]). As long as no vanishing events take place, 
that means as long as no small grains with topology class 2 < n < 5 do exist, switching events do likewise 
not occur, and ([7]) reduces to a system of uncoupled transport equations. This is in contrast to the 
gradient flow dynamics in which edge switching events occur independently of grain vanishing. 

2.3 Qualitative properties of the kinetic model 

We summarise the most important properties of the kinetic model. In order to simplify the presentation 
we argue by means of formal analysis but mention that all results will be proven rigorously within fJ31 

Decreasing number of grains Since equations ([7]) reflect a coarsening process, it is clear that the 
total number of grains N{f{t)) with 

oo 

should decrease in time. This is also satisfied for solutions to our model as grains with topology class 
n < 6 can shrink to area zero so that they are annihilated. More precisely, the evolution equations imply 

d ^ 

^^(/(i)) = ^ (n - 6) /„ (0, = ^ (n - 6) /„ (0, t) < . 

n>2 n=2 

Polyhedral formula In this section we motivate the choice of F in (l6|) . It is essential that our kinetic 
model reflects all the properties which are satisfied by a grain configuration which covers the complete 
area and where edges only meet in triple junctions. Hence we need to ensure that Euler's polyhedral 
formula is satisfied. For a finite network of grains Poincare's version of the polyhedral formula reads 

V + F-E = xig). 

with V, E, and F being the number of vertices, edges, and facets, respectively. Moreover, g is the genus 
of the surface and xid) — 2 — 2(7 the corresponding Euler characteristic. We can encode that grain 
boundaries only meet in triple junctions by setting V — 2/3E, and the polyhedral formula reduces to 

3F~E^3xi9)- 

In what follows we set xid) — because the network of grains is usually considered on a two-dimensional 
torus. In the kinetic model the normalised numbers of facets and edges are given by 

C30 oo 

F{f) = N{f) - E / E{f) = ^ E 

n>2 g n>2 q 
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respectively, and consequently we require each solution of ([7]) to satisfy 

Pifit)) = (9) 
for all times t > 0, where the polyhedral defect is given by 

oo 

P{f) = J2in-6) f .Ua)da. (10) 

n>2 g 

The main observation is that our choice of T guarantees (O. Indeed, due to ([5]) we find 

oo oo 

:P{f{t)) ^-Y.in-6f I djn (a, da + ^ (n - 6) 1 r (/ it)) (Jf)^ (a, t) da 

n>2 Q n>2 q 

oo 

= J2{n~ 6f fn (0, t) + T if (t)) Y^nj {Jf)^ (a, t) da, 



d 

It' 



i>2 



and a simple calculation shows 

E " iJf)n = PY.''^''- 1) /"-I + (/3 + 1) ^ n (n + 1) /„+i - nV„ - (/3 + 1) ^ 

n>2 Ti>3 n>2 n>2 n>3 

= 2(/3+I)/2-^n/„, 

n>2 

which implies 

thanks to (O. Hence, the polyhedral formula Q is satisfied for all t > if it is satisfied by the initial 
data. 

Conservation of area As a consequence of the polyhedral formula we obtain that the total covered 
area A{f{t)) with 

oo 

^(/) = E / «/"(«) da (11) 

n>2 Q 

is a conserved quantity. This follows from 



d 
dt 



OO OO 

Mfit)) = - E - 6) / adafn {a, t)da + r if (t)) JaY, {Jf)n («: da 

n>2 Q Q ""-^S 

oo 

= ^ (n - 6) / /„ (a, i) da = P{f{t)) = 0, 



where we used an integration by parts and the zero balance property 



2.4 Statement of the main result 



Our main result in this paper concerns the existence of mild solutions to the following initial and boundary 
value problem, provided that the initial data satisfy certain assumptions. 

Problem 2.1. For fixed <T < oo, and given initial data g = gnia) we seek mild solutions f = /^(a, t) 
to the following infinite system of coupled transport equations 

dtfn (a, t) + {n- 6) dafn (o, t) ^ T [f (t)) (Jf)^^ (a, t) 

with initial and boundary conditions 

/„(a,0) = 5„(a) for n>2, 

/„(0,t) =0 /or n>6, 

where a>0,te[0, T], n>2. 

A first assumption we have to make concerns the choice of /3. In what follows we always suppose that 
/3 is a fixed constant with < (3 < 2, where the upper bound is necessary in order to ensure that T]j{f) 
is non-negative for all / > 0. Further necessary assumptions regard the initial data. 

Assumption 2.2. Suppose that 

1. 5 is non-negative with 5n(0) = for all n > 6, 

oo 

2. J2 (l + "-)ll5«llo, 1 < with llffnllo, 1 = / {I + a)gn{a)da, 

n>2 ' 

3. g fulfils the polyhedral formula P{g) — 0. 

Moreover, suppose that g is sufficiently regular and decays sufficiently fast in n- and a-direction, 

According to the discussion in f}2l the first three assumptions on the initial data appear very naturally. 
Our regularity and decay assumptions, however, are needed for technical reasons and can probably be 
weakened at the price of more analytical effort. The precise statement of these assumptions appears below, 
but we mention that we mainly assume all functions gn to be equi-continuous and to decay exponentially 
with respect to a and n. 

Our main result can be summarised as follows. 

Theorem 2.3. For any initial data that satisfy Assumption \2.S\ there exists a unique mild solution f to 
Problem p.ip for all <t < oo. Moreover, this solution conserves the area with non-increasing number 
of grains, and all states f[t) satisfy Assumvtion \2/A 

The details of the proof are presented within |J3] and rely on the following key ideas. In H3.ll we 
introduce an approximate system of coupled transport equations by neglecting all topological classes 
with n> uq. 

In order to construct mild solutions for the approximate system we solve the transport equations 
explicitly, and apply Duhamel's Principle (or Variation of Constants) to the collision operator. This is 
discussed in ij3.2l 

Freezing the coupling weight F in the approximate equations we can derive a comparison principle 
that implies both, the non-negativity of solutions and the existence of an appropriate super-solution. 
This will be done in p. 31 

In §3.4l we establish the short time existence and uniqueness for admissible solutions to the approximate 
system. To this end we propose a suitable iteration scheme, and make use of Banach's Contraction 
Principle. 
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The decay behaviour with respect to a and n is investigated in H3.5I We derive several tightness 
estimates, and as a consequence we obtain long-time existence for the approximate system together with 
estimates that are uniform in tiq. 

In i j3.6l we show the existence of mild solutions to Problem 12. II bv passing to the limit no ^ oo. The 
Arzela-Ascoli theorem provides that solutions to the approximate problem converge to some reasonable 
limit, and the tightness estimates ensure that this limit provides an admissible solution. Finally, we 
sketch how both, the uniqueness of solutions and the continuous dependence on the initial data, can be 
obtained. 

2.5 Further considerations 

Here we point out some open questions and directions for future research. 

Stationary solutions Nontrivial stationary solutions to (O, i.e. ^ 0, are characterised by (a) = 
for n ^ 6 and < a < oo (cf. [Hen07| . §6.1), but the component function /| (a) can be anything. This 
is a further difference to the gradient flow dynamics as there not all networks of hexagons are stationary. 
However, such nontrivial stationary solutions to the kinetic model are expected to be unstable for the 
following reasons: Slightly perturbed data lead to a positive F (/ (t)) for some times t and are therefore 
affected by the coupling operator {Jf)^ {a,t). This leads to a decrease of the total number of grains. We 
therefore expect stationary solutions to be unstable, and thus we suspect that / ^ weakly for t oo. 

Self— similar scaling For a coarsening process as considered here, one usually expects to find self- 
similarity under dynamic scaling. The natural rescaling ipn (^, t) = i^/„ (a, t), ^ = a/t, yields the following 
equation for self-similar solutions 

(n - 6 - ^cv^n - F (^) ( J(^)„ + 2(^„ , n>2, ^>0. (12) 

The natural boundary conditions are ipn (0) = for n > 6, so that the solution depends on the values of 
(/32 (0) , . . . , (/?5 (0). Note that the coupling weight still depends on the complete solution. A starting point 
for future analysis is the following observation. We can integrate ([T2l) with respect to ^ to obtain 

(6 - n) ^„ (0) = F (^) ( J0)„ + 0„ , n>2. 

This is a two-point iteration scheme for = J (C) d^- 

Lewis' law A natural question concerning grain growth is to ask whether there are correlations between 
the topological class and the area of a grain. Lewis [Lew43| observed a linear relationship examining 
cellular structures arising in biology, and Rivier and Lissowski |RL82) derived Lewis ' law by maximum 
entropy arguments applied to cell distributions. In common with Flyvbjerg [Fly93| this so-called Lewis' 
law reads 

(0„ = &("-6)+c (13) 

in our model. Here = / 5(p„ (5) d^/ / {£_) d£_ denotes the mean grain size in the topological class 
n. However, it is unclear if this phenomenological law is really applicable for grain growth. Formal 
computations suggest in our model that (fT^ is valid asymptotically for large n with b = 1/ (F + 1) and 
c = 6 ((2/3+ 1) - 6F) (cf. |Hen07| . §6.4). Similar results are achieved by Flyvbjerg |Fly93| . 

3 Proof of the main result 

3.1 The approximate system 

The approximate system for ([7]) results from the original equations by neglecting all number densities 
belonging to topological classes with n > uq. More precisely, we choose the parameter uq with 6 < no < 
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cxD, and modify the coupling operator accordingly. The approximate coupling operator J splits into its 
gain and loss part, that is 

J = J+-J_, (14) 

which now are given by 

(J+/)2 = 3(/3+l)/3, 

(J+/)„ = (/3 + l)(n+l)/„+i+/3(n -!)/„_! for 2 < n < np, (15) 
(J+/)„^ =/3 (no -1) /„„-!, 

and 

(J_/)2 = 2/3/2, 

(J_/)„ = (2/3+1) n/„ for2<n<no, (16) 
(J_/)„^-(/3+l)no/„„. 

Notice that for n < ng the term ( J/)„ is defined as in the original model. For the sake of consistency we 
must moreover adapt the formula for F. In what follows we use the approximate coupling weight 

r(/)-rA.(/)/Fz,(/), F^(/) = ^(n-6)'/„(0), TD{f) = -Y. f n{Jf)^{a)Aa. (17) 

Analogously to [J5]we define the total area A(f), the number of grains N{f), and the polyhedral defect 
PU) by 

no «o °? no °? 

^(/) = E / /"(«) A{f) = J2 afnia) da, P(/) = ^ (n - 6) / /„(a) da, 
and for convenience we introduce in addition 

oo oo oo 

M{f) = Y,^ J /n(a)da, i?(/) - 2(/3 + 1) j /2(a) da -no/3 j /„„(«) da. 

n=2 Q 

Remark 3.1. 1. Definitions (fT5|) and (fT6|l imply 

no no no 

E("^/)n = 0' E"(^/)n = 2(/3 + l)/2-no/3/„„-^n/„, (18) 

n=2 n=2 n=2 

and we infer that 

AiJf) = NiJf) = 0, P( J/) = Af ( J/) = -F,,(/). (19) 
2. The polyhedral formula P{f) = implies 

rD{f)^-R{f) + 6N{f), 

and if / is in addition non-negative we have, thanks to < /3 < 2, 

rD(/) > (6-2(/3+l))7V(/) >0. (20) 
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The approximate problem we aim to solve can now be stated as follows. 

Problem 3.2. For fixed uq > 6 and < T < oo, and given initial data g = gnia) we. seek (mild) solutions 
f = fn{a,t) to 

dtfn (a, t) + {n- 6) dafn (a, t) = r (/ (t)) iJf)^ (a, t) (21) 

with initial and boundary conditions 

/„(a,0) = g„(a) for 2 < n < Uq, , . 

fn{0,t) =0 for 7<n<no, ^ ' 

where a>0,t€z[0,T],2<n< uq, and J and T are defined as in (jl5p - (|17p . 
3.2 Transport equations and notion of mild solutions 

We start with some basic facts about solutions to transport equations in the upper-right space-time 
quadrant a > and < > 0. For this reason let us consider the following system of transport equations. 

Problem 3.3. Let a > 0, 2 < n < uq and t G [0, T]. For fixed initial data g = gn{o) and given right 
hand side h = hn[a,t) we seek (mild) solutions f = /„(a,i) to 

dtfn + {n- 6)dafn = K (23) 

with initial and boundary conditions as in the approximate problem, see (j22p . 

The homogeneous problem with h = Q can be solved explicitly by the method of characteristics. This 
means, the general solution to the homogeneous problem is given by f[t) — T(t)g, where the group of 
transport operators is defined by 

with 

(T„_6(i)5n)(a) = 9n[a - {n - 6)t) 
for all n < 6, whereas n > 6 corresponds to 



gn{o, — (n — 6)t) for a > (n — 6)t, 
for a < (n - 6)t. 



Recall that n < 6 implies the transport velocity n — 6 to be non-positive, so there is no contribution from 
the boundary in this case. 

For non- vanishing right hand side h the solutions can be constructed by means of Duhamel's Principle. 
More precisely, for given h = hn{a,t) the unique mild solution to Problem 13. 31 is given by 

t 

fit) = T{t)g + J^it- ds. (24) 



Moreover, the mild solution depends continuously on the data via 

t t 
||/(t)IL < ||5lL + J IIMs)ILds, ||/(i)llo,i <C^no.*(ll5llo,i+ / ||Ms)llo,ids), (25) 



where \\g\\o,i ■= T,7=2 1^ (1 + a)l5«(a)|da < oo. 

In the sequel we make use of the following regularity results for mild solutions, that can be derived directly 
from the representation formula given in Ij24p . 
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Definition 3.4. 1. A state g — g„(a) is called 

(a) C'^ -regular, if g is continuous (w.r.t. to a) and — gn{0) for all 7 < n < no, 

(b) L^-regular, if H^Hq ^ < oo. 

2. The right hand side h = hn{a, t) is called 

(a) C'^ -regular, if h is continuous (w.r.t. to (a,t)), 

(b) L^-regular, if ||^(^)|lo ^ dt < oo. 

Moreover, we say the solution / = /n(a, t) is regular, if / is a regular right hand side and f{t) is a regular 
state for all t e [0, T], where 'regular' means either C*^- or L^-regular. 

Lemma 3.5. The mild solution f to Problem \3.3\ given in (j24p has the following properties. 

1. If the data are -regular, then f is -regular. In particular we have fn{0,t) — for all t and 
7 < n < uq. 

2. If the data are -regular, then f is -regular. 

Finally, we collect some properties of mild solutions to be used in ^'S.4[ 

Remark 3.6. If the data g and h are compactly supported in {a : < a < qq}, then the mild solution 
f{t) is compactly supported in {a : < a < oq + (no — 6)^} 

Lemma 3.7. Suppose that the data {g, h) for Problem \S.S\ are C*^- and -regular. Then the mild solution 
f satisfies 

* 5 * 

Nifit))^Nig)+ J iV(/^(s))ds + ^(n-6) J /„(0,s)ds, 



t t 

P{f{t)) = P{g) + j P{h{s)) ds + J Toifis)) ds, (26) 





t t 



A{f{t))=Aig) + J Aih{s))ds + J P{fis))ds. 



In particular, N{f{t)), P{f{t)), and A{f{t)) are continuously differentiahle with respect to t. 

Proof. For classical solutions the differential counterparts of all assertion follow directly from the dif- 
ferential equation (I23|) . For mild solutions we approximate with classical solutions. In fact, we can 
approximate the data {g, h) by some regular data [g, h) that satisfy the boundary condition and are both 
differentiahle and L^-regular. The corresponding solution / to Problem 13.31 is then a classical solution 
and satisfies ((2^ . Finally, / — / is a mild solution to Problem 13.31 with data {g — g,h — h), and using 
(US]) we infer that / satisfies □ 



3.3 Auxiliary problem with prescribed coupling weight F 

In this section we consider an auxiliary problem that results from the approximate svstem [3T2l bv prescrib- 
ing the coupling weight F as a function of time, and study both the existence and qualitative properties 
of solutions in the space of bounded and continuous functions. In particular, we will derive comparison 
results for such solutions. 
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Below in ij3.4l we apply our results in the context of the approximate problem, and show that each solu- 
tions to the approximate problem must be non-negative and bounded from above by some appropriately 
chosen super-solution. 

The auxiliary problem can be stated as follows. 

Problem 3.8. Let < T < oo be arbitrary, T g C([0, T]) be some non-negative weight function, and g 
be some initial data. Then we seek (mild) solutions f to 

dtfn + {n-6) dafn = r(t)(J/)„ (27) 

with initial and boundary conditions (j22p . 

The state space and solution space for Problem 13. 81 are given by 

= BC{[0, oo); Xaux) and J^aux = C([0, T]; Faux), 

respectively, where Xaux = K"°^^ denotes the set of all tuples {x2, ...,Xno)- 

The key idea for arriving at comparison results for the auxiliary problem is to split the exchange operator 
into its loss and gain part according to (fH|) . More precisely, we regard each solution to the approximate 
problem ((27l) as mild solution to 

dtfn + (n - 6) dafn + r(i) J-/(i) = Ht) , (28) 
where the right hand side is given by 

h{t) = r{t)j+f{t). (29) 

Doing so we benefit from the following two observations. The loss operator J_ is diagonal and can 
hence easily be incorporated into the homogeneous problem. The resulting solution operators take care 
of the transport as before, but describe additionally the exponential relaxation of /„ to along each 
characteristic line with local rate proportional to T{t). In particular, these solution operators preserve 
the non-negativity of the initial data for all times. On the other side, the gain operator J_|_ preserves 
the cone of non-negative functions, so that the non-negativity of solutions turns out to be a direct 
consequence of Duhamel's principle with respect to the gain operator. 



Exponential relaxation along characteristics 

The mild solutions to the homogeneous problem corresponding to ([25]) are described by a two-parameter 
family of linear operators defined by 

t 

Tf{s,t) : Yaux Kiux, {Tf{s,t)f)^ := exp ^-m„ J T{s)dSj Tn-&{i - s)fn, (30) 

s 

where s and t are two times with < s < i < T, and w„ is the relaxation constant appearing in the loss 
operator. More precisely, according to (|15p and ([1^1 we have 

U2 — 2(3, u„„ = (/? + 1)^0, and Un — {2(3+l)n for 2 < n < uq. (31) 

With ((Sn)) . each mild solution to satisfies 

t 

fit)^Tf{0,t)g + J Tfis,t)his)ds, 



with h as in (|29p . As outlined above, the introduction of 7^~(s,i) is motivated by technical reasons but 
does not change the auxiliary problem. More precisely, we find the following equivalence. 
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Remark 3.9. A mild solution / e 3^aux to Problem 13.81 that is 



t 



fit) = T{t)g + Tit- s)Tis)iJfis)) ds, 



also satisfies 







t 

/(t) = Tp-(0,t).9 + J Tfis,t)Tis)iJ+fis))ds, 





and vice versa. 

Super solution to the auxiliary problem 

In order to establish an existence result for mild solutions we start by identifying a suitable super-solution 
(/) to Problem 13.81 that does not depend on t and a, and provides suitable a priori bounds, see Lemma 

EH 

Remark 3.10. Let G Xng be defined by 



4>nia) 



Then, is a solution to J(j) = 0. Moreover, / with /„(a,t) = 4>n for all t > 0, a > and 2 < n < uq 
provides a solution to the differential equations (ITfl) . 



Proof. For n > 3 we find 

,„:=(/3 + l)n0„-/3(n-l)0„.,= ^^j^^j - j =0, 

and this implies (J0)„ = jn+i — jn = for all 2 < n < ng. Moreover, we have (J0)2 = is = and 

(•/'^)„o = -J"o =0. □ 

We next show that a suitably chosen multiple of / provides an upper bound for each solution to the 
approximate problem. For this reason we introduce a special norm for tuples ix2, ...,a;„o) G Xaux by 

\x\a := sup (32) 

2<r!<no + l 'Pn 

This norm can easily be extended to states / € Yaux via ||/||^ :— supj,>g |/(a)|^, and writing 

11/11/3= BUp L/J„,^, L/J„,/3:=SUp^^ (33) 

2<ri<no + l a>0 <Pn 

we infer that 

ll/„IL<ll/IIA^^/3||/||^. 

Notice that the norms || • ||^ and || ■ are equivalent on Yaux since tiq < oo. However, we prefer working 
with II ■ |L because only this norm gives rise to estimates that are independent of no. 
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Remark 3.11. For all / 6 Yaux and 2 < n < no we have 

||(J+/)JL<«n</>„||/||0 

with Un as in (PT|) . 

Proof. It is sufBcient to consider ||/||^ 1. For 2 < n < Uq we find 
||(J+/)JL<(/3+l)(n + l)||/„+i|U+/3(n-l)||/„_i||, 



Finally, the estimates 



and 



(•^+/)nolloo = /^("O - l)||/no-l|loo < ( ^-pyj ^ {jTl) ) ^ ^^0<t>n 



||(J+/)2|L = (1 +/?)3||/3lL < ^f^V - = «2</'2 



complete the proof. □ 
Existence of solutions to the auxiliary problem 

Our existence and uniqueness results for the auxiliary problem are based on Banach's contraction principle 
applied to the iteration operator 

t 

laux : 3^aux ^ 3^aux, (2'aux/)(t) = {t)g + j Tf{s, t)T [s) {J + f [s)) ds. 



We start with deriving some a priori estimates for 2aux- 
Lemma 3.12. For each f e J^aux we have 

sup ||Zaux/(i)ll/3 < max{||.g|L, sup \\f{t)\\A. 

0<t<T 0<t<T 
t 

Proof. Let u„ be as in (PT|) and G{t) = /r(s)ds. Definition ([HU]) and Remark l3 . 1 II imply 



||(Xaux/)„WIL <'^nfexp(-u„G(0)||3||^+ / exp(-u„G(0+^i„G(s)Kr(s)||/(s)||^ds 
and we conclude 

L(2'aux/)(0J„,/3 < exp (-u„G(i)) ^||g||^ + ^(t) y exp (u„G(s))u„r(s) ds 

= exp(-u„G(i))(||5||^ + /i(t)(expKG(t)) - 1)), (34) 
with = supo<s<t ||/(s)||0, which yields the desired estimate. □ 
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We mention that all results obtained so far in this section immediately apply to the case no = oo, 
provided that the initial data satisfy ||5||^ < oo. 

Corollary 3.13. For arbitrary < T < oo and any initial data g £ i^aux there exists a unique mild 
solution f e 3^aux to Problem \3.8\ which satisfies ||/(t)||^ < llsll/j /o*" t G [0, T]. Moreover, this 
solutions is 

1. non-negative provided that the initial data are non-negative, and 

2. C^-regular (L^-regular) if the initial data are C^-regular (L^ -regular) . 

Proof. The existence and uniqueness of a mild solution is a direct consequence of Banach's contraction 
principle and Lemma 13.121 (the Lipschitz constant of Xaux with respect to || ■ ||^ can be read off from p4p 
with 5 = 0, and is given by 1 — exp (— u„(,G(t))). Moreover, the remaining assertions are satisfied as Xaux 
respects both the non-negativity and -regularity. □ 

3.4 Solutions to the approximate system 

In this section we are going to establish local existence and uniqueness results for the approximate 
problem. To this end we consider the solution space J^app — ^^([0, T]\ Yapp), where the state space given 

by 

^app = ^aux n {/ = /„(a) : ||/„||q ^ < oo for all 2 < n, < no}. 

Of course, the set of all admissible states is only a proper subset of llpp as any reasonable solution to 
the approximate problem must be non-negative and area conserving, and has to satisfy the boundary 
conditions and the polyhedral formula. 

Definition 3.14. The state / G Yapp is called admissible, if it is non-negative with A{f) > 0, and satisfies 
the polyhedral formula P{f) = as well as the boundary conditions /n(0) = for n > 7. 

Iteration scheme for the approximate system 

We construct solutions to the approximate system as fixed points to the following iteration operator 
^'app : J^app J^app- For givcu / e 3^app let / = lapp/ be the mild solution to 

dtfn + {n-6)dafn^\r{f)\{Jf)^ (35) 
with initial and boundary conditions (|22|) . 

Lemma 3.15. For all initial data g 6 Yapp there exist a time T > and a constant Co (both depending 
on no and g ) such that the operator Xapp has the following properties. 

1. For given f £ ^aux the function f — lapp/ G 3^aux is a mild solution to (1351) that satisfies the initial 
and boundary conditions (|22p . 

2. For -regular data g and f the function J^ppf is -regular. 

3. The set y,pp = BC{[0, T]; Y,pp) with 

?app = {/ e Fapp : r,,(/) > C^-\ ll/ll^ < Co} . (36) 

is invariant under the action o/Iapp, and the operator f |r(/)|J/ is Lipschitz- continuous on 
Y 

app- 
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Proof. At first let Cq and T be arbitrary but fixed, and suppose tliat f{t) e Yapp for all t G [0, T]. For 
the remainder of this proof C always denotes a constant depending only on uq and /3, but the value of C 
may change from line to line. From (|35p and (|25p we conclude that 

t t 2 



and due to / G Yapp we have 

ll/»lloo<ll5lloo+^C'o'r. (37) 
Moreover, since F^i is linear in / we find 

t 

TD{.f{t))~TD{mg)\< J r^l^j^j)' \^D{T{t - s)Jf{s))\ ds < CC.^T. (38) 



In order to prove the invariance of Yapp under the action of 2app, we choose Co sufficiently large such 
that 

2||<7lL<C^o, Fz3(g)>3Co-\ 
In addition, we choose T sufficiently small with 

ri3(r(t)5)>2Co-^foraUO<t<T, CCq'T < C C^^ T < 

and thanks to (|37p and we thus find / £ 3^app- Finally, the claimed regularity results are provided 
by Lemma 13.51 ^-nd the proof of the Lipschitz-continuity is straightforward. □ 



Local existence and uniqueness 

We prove the existence of admissible solutions to the approximate system by combining the following 
arguments: (1) Banach's contraction principle implies the existence of a unique fixed point / for Zapp- 
(2) Our results concerning the auxiliary problem provide the desired non-negativity. (3) The polyhedral 
formula is a consequence of the equation itself, and implies the conservation of area. 

Lemma 3.16. Let T and Cq he as in Lemma \3.15\ and suppose that the initial data g are admissible. 
Then there exists a unique mild solution f e 3^app to the approximate system on the time interval [0, T]. 
Moreover, 

1. all states f{t) are admissible (in the sense of Definition \3. 14\ l and fulfil \\f{t)\\p < HffH^, 

2. f conserves the total area and the number of grains is non-increasing. 

Proof. Banach's contraction principle provides the existence of a fixed point / G 3^app with I^ppf = /. 
Then we consider the auxiliary problem with prescribed coupling weight |F(/(i))|, and both the positivity 
and boundedness of / are consequences of Corollarv 13.131 In particular, we find F(/(i)) > and this 
shows that / is indeed a solution to the approximate system. In order to prove the polyhedral formula 
and the conservation of area we use basically the same computations as in ^ More precisely. Lemma 
13.71 implies 

1 no 

-JlPifW) - E (" - 6)'/(0,t) = T{f{t))P{Jf{t)l 

n=2 
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and with (fT7)) and we infer that 



Ap(;(,)).r.(/(t))-M^r.(/W) = o. 

This imphes P{f{t)) — for aU t due to P(g) = 0. Similarly, we find 

±A{f{t)) = P{f{t))=Q. 



Finally, exploiting Lemma |3 . 71 once again, gives 

-N{f{t)) = E - 6)/n(0, t) + V{f{t))N{Jf{t)) < 



dt 



where we used N{Jf) = and / > 0. □ 

Remark 3.17. The proofs of Lemma [3.151 and Lemma [3.161 imply that the local solution does exist 
as long as T]j(f{t)) can be bounded from below by a positive constant da. In fact, for non- negative 
/ both and rAr(/(<)) are bounded by the super-solution f^{a,t) — ||ff||^</'ra, and thus the 

constant Cq from (j36p can be chosen to depend on do Sind g only. Consequently, each lower estimate for 
info<t<T rD(/(i)) will imply the existence and uniqueness of solutions on the interval [0, T]. In what 
follows we derive several variants of such estimates by bounding N(f(t)) from below, respectively, which 
is sufficient in view of (PO]) . 

The first lower bound for the number of grains is an elementary consequence of the von Neumann- 
MuUins law. The main observation is that grains can grow only with speed hq — 6, so that for compactly 
supported initial data the number of grains can be bounded from below by the total area. More precisely, 
if the initial data g are compactly supported in {a : < a < uq} then we have 

Nifit)) > /Ij-'^ (39) 

This bound is, however, not optimal because it holds only for compactly supported initial data and 
depends strongly on uq. For this reason we finally rely on refined estimates to be derived in Lemma [3.25l 



3.5 Tightness estimates and global existence 

We introduce the notion of 'quasi-complement', that is the number of grains outside a bounding frame. By 
enlarging this frame in time, we show that no mass runs off at infinity in finite time, and as a consequence 
we establish long-time existence for the approximate system. Moreover, these tightness estimates also 
apply to the original problem with uq — oo provided that the initial data decay sufficiently fast with 
respect to both a and n. 

Within this subsection we use the following notations. Let / be a fixed non-negative solution to the 
approximate system (|3.2p . defined for < < < T < oo, and set 

7 = iog(n- 1//3), Nit)^N{f{t)), r(<)-r(/(i)), r(t) = sup r(s). 

0<s<t 

Notice that the definition of 7 implies Pn<j)n — exp (— ^7) with <j) = {4'n)2<n<no being the super-solution 
from Remark 13.101 

Given G N and a G [0, 00) we define the first quasi-complement by 

N^{t,a,v)= fnia,t)da + J2 fnia,t)da 
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as the non-essential part of N (t). In order to simplify the notation we allow for = by setting /i = 
and (J/)^ = 0, so that N{f{t)) = N^{t,0,0). 

In what follows we write {t, a) := {t, a, [a]) if the third argument of is given by rounding 
down the second one, and refer to 

>i oo ^, CXD 

no ^ no ^ 

M^(i,a)= J2 n U{a,t)da + J2 afn{a,t)da (40) 

n=LaJ+l { n=2^ 

as the second quasi-complement. 

Remark 3.18. For aU t > and a > we have N^{t, a) < M^{t, a). 

Proof. The desired estimate follows immediately from the definitions provided that a > 1 . For < a < 1 
we find 



oo 

N^{t, a) = N{t) <^n f„{a, t) da < M^{t, , 

n "J 



due to [aj =0. □ 
Estimates for N-^ 

In order to derive suitable estimates for the first quasi-complement we choose v = and a as 
functions of time. More precisely, a should grow at least linearly in v to compensate for the transport in 
a along characteristic lines, and v should grow exponentially to control the diffusion in n. 

Lemma 3.19. Let /z, a : [0, T] [0, oo] be two smooth functions with 

fi{t) > and a{t) > max{ l^l{t)\ - 6, 0}. 

Then, we have 



{t,a{t) ,l^I{t)\) < {0,aiO) ,[^iiO)\) + \\g\\^T{t) J exp (-7L/i (s)J) a (s) ds (41) 



for all t e [0, T] . 

Proof. We define v {t) = [/x (i)J as the integer part of ^jL{t) and denote the jump of v [t) by \y\{t) = 
v{t+) — v{t—). Obviously, we can restrict ourselves to v{t) < uq, and for simplicity we consider classical 
solutions. Our results then can be generalised to mild solutions by approximation arguments. At first we 
study the case {vjit) = 0. Differentiating N-^ {t, a{t) ,v {t)) and using the evolution equation (PT|) yields 

d 



N^it,a,,y) = - J2 {n-6)fnia,t)+T{t) ^ / ( J/)„ (a, i) da 

n— f+1 n— Q 

no no °9 

a/n(a,i) / (n-6)5a/n(a,t)da (42) 

— ,,11 ^ — n ^ 



no 



^(J/)„(a,t)da-^a/„(a,i), 

„ n o 



n=2 n=2 
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dt' 



where a and v are shorthand for a(t) and ^{t), and due to psp this impHes 

-7V^ (i, a, z.) = ^ (n - 6 - «) /„ (a, t) + T{t) ^ / ( J/)„ (a, i) da 

n— 2 n— Q 

< / (J/)„(a,t)da 

For 1/ = or = 1 we find ^N-^{t, a) < due to X;"l2 (-^Z)™ = 0' "^^^^^ for > 2 the identities ^ 
and (fTC|l give 

"0 

< (a,i) < 1 15 1 1/3/3 i^'/'n = ||5||^exp(-7:/), 
where we used that ||/(i)||^ < \\g\\p, see Corollary 13. 131 Therefore, 

^AT^ (i, a{t), j.(t)) < r(i) I \g\\p exp (-7^. (t)) a (t) (43) 

holds for all v. In the case = +1 we find 

IN^ {t, a, La^J ) 1 = (t+, a(t+), [mJ (^-) + I) - {t-,a{t_), {t.}) 
«(*) 

f[Kt-}\ da < 

as an additional part in the r.h.s. of so (I43p still holds in a distributional sense. Finally, (|4ip follows 
by integrating (jiSj) . □ 



Next we derive a decay result for iV^ {t, a) with respect to the variable a which does not depend on 
riQ but only on the initial data g. These decay estimates turn out to be the crucial ingredient for both, 
the global existence proof for solutions and the passage to the limit no — > 00. 

Corollary 3.20. The estimate 

{t, a) < (0, a exp (-t)) + C\ I5I |^ T{t) exp (-7 a exp {-t)) 
holds for all t > and a > 0, where C depends only on [3. 

Proof. Let > be arbitrary and consider a{t) — ^{t) — aoexp(t), and let i' (t) = [/i(t)J. Then, 
a{t) ~ ii{t) > h'{t) — 6, and Lemma 13.191 implies 

t 

{t,aoexp{t))<N^ (0, ao) + \\g\\f^T{t) exp (7) J exp {-jfi (s)) a {s) ds, 



where we used that exp (— 7[/i (s)J) < exp (—7/1(5) +7). Replacing s by a = a{s) = Q;oexp(s) we find 

ao oxp (t) 

7V^(t,aoexp(i)) <7V^(0,ao) + ||g||^r(t) exp(7) J exp {-j a) da 

ao 

< (0, ao) + \\g\\^ T{t) exp (-7 ao) . 
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Since this identity holds for arbitrary t and ao, we can choose ao — aexp {—t), which gives the desired 
result. □ 

Remark 3.21. The proofs of Lemma lS.lQI and Corollarv l3.201 and hence all estimates derived below, can 
be easily generalised to the original problem with rto = oOj provided that WgW^ < oo, i.e., if for n —^ oo 
the term ||5n|loo decays at least as fast as n"^ cxp (—771). 

Estimates for 

Here we exploit the properties of the A''-'- and that the decay behaviour of Af-^ with respect to a is 
controlled by t and the first quasi-complement of the initial data. 

Lemma 3.22. We have 

it, a)<C exp (t) (l + | |g| |^ + | \g\\^ r(i)) A/'o^(a cxp (-t)) 
for all t > 0, a > 0, where A/'q" depends only on the initial quasi-complement via 

oo 

Af^{a) = (a + l) exp{--fa) + aN^ {0,a) + ^ iV^ (0, a) da, (44) 

a 

and C depends only on (3. 
Proof. Notice that 

„ oo ^, oo oo ^, oo 

no ^ ^10 /. /. f> 

/ a/„ (a, i) da = ^ / / /„ (s, i) ds da + a ^ //„(a,i)da (45) 
via an integration by parts. CoroUarv 13 . 201 vields 

oo 

/ /n(a,t)da< A^^(0,aexp(-<))+ C||g||^r(i) cxp(-7acxp(-i)), (46) 
and this bounds the second term on the r.h.s of (|45|) . Moreover, integrating (|46|) we find 

oo oo oo 

51/ / /n (s,t) dsda < + exp(i) / iV^ (0, a) da 

a Qcxp(-t) (47) 

C||5||;3r(i)exp(t) 



7 

and combining with ([T7|) we conclude that 



■ exp (—70; exp (— i)). 



^ / a/„(a,t)da<Cexp(t)(l + ||5||^ + ||5||^r(<))AAo^(acxp(-0). 

n=2^ 

To derive the remaining estimate we split the integration with respect to a at a = n, and obtain 

00 n 00 

no r- no „ no „ 

n fn{a,t)da< ^ n fn{a,t)da+ ^ I afn{a,t)da 

„=[qJ+1 p n=[Qj+l Q n=[Qj+l„ 

no "0 °? 

n fn{a,t)da + 'Y afn (a, i) da. 



< 

4aJ+l ^ n=2; 
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We have already bounded the second term on the r.h.s., while the first one can be estimated by using 
fnia,t) < ||5||^0„ = ||5||^exp(-77i)/(n/3). This gives 



" y /n (a,i) da < -^^^^ ^ nexp (-771) < C ||g||^ J aexp{-ja)da 

n=la]+l Q n=[a\+l ^ 

00 

(aexp {-t)) 

a cxp {—t) 

and the proof is complete. □ 
Positivity of numbers of grains and global existence 

The initial data g are called rapidly decreasing if there exists two constants do and Dq such that 

M^(0, a) < Do exp (-do a) < 00 (48) 

for all a > 0. Notice that (PS)) implies that both N^(Q,a) and Af(^{a) decay exponentially with respect 
to a. This follows from Remark 13.181 and Definition 



Remark 3.23. Suppose that the initial data decay exponentially with respect to both a and n, this 
means there constants do and Do such that gn{ci) < Doexp (—do{a + n)). Then there exists a suitable 
choice of do and Dq such that is satisfied. 



For rapidly decreasing initial data we can estimate M'^{t, a) for arbitrary t and a by means of Lemma 

EM 

Corollary 3.24. Suppose that the initial data g are rapidly decreasing. Then, for each t > there exists 
constants dt and Dt such that 

M^{t,a) < Dt{l+T{t)) exp{~dta) 

holds for all a > 0. These constants depend on (3 and g, hut not on T{t) or no- 

Recall that the lower bound for N{t) from ((5^ is not optimal, as it is restricted to compactly supported 
initial data and depends on no. Here we derive a better result by exploiting the tightness estimates. 

Lemma 3.25. For rapidly decreasing initial data g and each i > there exists a constant Ct > such 
that 

N(t) > Cu r(t) < i/Ct. 

In particular, this constant is independent of no . 

Proof. Within this proof let Ct denote an arbitrary constant that depends only on i, /3, and g. According 
to Remark |3. II and the monotonicity of N{t) we have 

r(i)= sup L s = sup -rr^ < -rr-r- 49 

o<s<t o<s<t N{s) N[t) 

Moreover, the conservation of area implies 

^ a 00 

no ^ /. 

A{g) = A{f{t))^J2{j afn{a,t)da + J a fn{a,t) da^ < aN{t) + M^{t, a) 
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for all a > 0, so we infer from Corollary 13 . 241 that 

aNit)> {A{g)- Dt{l+Tit))cxp{-dta)) > (^(5) - (l + Ar(i)"') exp (-dta)) . 
Now we choose a = at such that N{t)at = ^A{g) — ct, that means 

at = Ct{l + ln{l + N{t)-^^y 

and hence we find 

Ct (l + In (l + N{t)-^^^N{t) = 1. 

This estimate yields the existence of a constant Ct with N{t) > Ct, and thanks to (jH]) this provides also 
an corresponding upper bound for T{t). □ 

Corollary 3.26. Suppose that the initial data g are rapidly decreasing. Then the solution from Lemma 
\3.16\ does exist for all times t>Q. 

Proof. Lemma [3.251 provides a priori estimate for both N{t) and r(<) for arbitrary large t. Therefore, the 
local solution to the approximate system exists for all times, compare Remark 13 .171 □ 

3.6 Passage to the limit no — > 00 

In this section we consider a fixed final time < T < 00 and pass to the limit riQ 00. Consequently, 
from now on (i) the formulas for Jf and r(/) = ^N{f)/^D{f) refer to (jl]) and the definitions for the 
original problem with hq — 00, and (ii) the number of grains N{f), the area A{f), and the polyhedral 
defect P{f) are given by ([8|), (fTTj) . and (fTO]) . respectively. 

Existence of solutions 

The solution space 3^ for Problem 12. H is given y = C([0, T);Y) for some < T < 00 with state space 

Y:=[feBC{[0,^);X) : 1 1/| |^ + 1 1/| |„,o„. < ^} , 

where X = R^"--^^ abbreviates the space of real-valued series (x2, ...,a;„, ...). The norm || ■ ||^ is defined 
analogously to ((5^ and ([55)) (with no = 00), and || ■ \ \^^^ reads 

00 

ll/ILom:=E / (" + «)l/n(«)|da. 

In what follows we call a state f & Y admissible, if / is non-negative with A{f) > 0, and satisfies the 
polyhedral formula P{f) = and the boundary conditions fn{a) = for n > 6. 

Our strategy for constructing mild solutions to Problem 12.11 is rather straightforward. For fixed < 
T < 00 and given admissible initial data g G Y we consider a sequence of corresponding solutions to the 
approximate problem with increasing parameter no, and aim to show the existence of an reasonable limit 
in y. More precisely, within this section we consider the functions 

gM ^Y, /(""^ej^, r("°) e C([o, T]), 7v("«) e c([o, r]), m^^("«) e c([o, t] x [o, t]). 

which are defined for uq > 6 as follows. 
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1. For all a > and < t < T we have 



_ / foi" 2 < n < no, Hn„)(^ f fn{a,t) for 2 < n < no, 

^" ^ ' \0 for 71 > no, ^ ' ' \0 forn>no. 



2. For each n-o the component functions f^°\ /n""'' ^-re the unique mild solution for the approximate 
Problem 13.21 with initial data given by g^°\ 

3. r("") and iV("") arc the corresponding coupling weight and number of grains, respectively. 

4. M-*-^'"") is the second quasi-complement from 

Recall that Lemma [3.151 combined with Corollarv l3.26l provide that all these functions are well defined, 
and that each F^"") is positive and bounded. 

It is natural to suppose the initial data to be admissible, but we need a bit more for our subsequent 
analysis. In order to establish compactness in y and to control the tail behaviour (w.r.t. to n and a) we 
must assume that the initial data are moreover regular. 

Definition 3.27. A state / e F is called regular if the following statements are satisfied. 

1. / is equi-continuous with respect to || • ||^, i.e., for each e > there exists 5 = 5{e) such that 

\\Stf-f\\p<e for all < <5 < J, 

where is the left-shift operator defined by S'^ f{a) = f{a — 6) for a > S and S'g f{a) = else. 

2. / is rapidly decreasing w.r.t n and a, that means there exist two constants d and D such that 

M^{f,a) < Dexp{-da) 

for all a > 0, where 



oo oo 

M-'-(/, a):=^ / af,i{a)da+ ^ n /„(a)da 



n>2„ n>LQj+l 

is defined in line with (HO]) . 

We proceed with some remarks concerning the regularity of the initial data g. {i) The first condition 
is apparently satisfied ii g is differentiable w.r.t a with ||9a(?||^ < oo. (ii) Concerning the second condition 
recall that \\g\\g < oo already implies an exponential decay w.r.t. n. (in) Compactly supported initial 
data with gn (a) = for large n and a are rapidly decreasing. 

Assumption 3.28. From now on we suppose the initial data g e F to be admissible and regular. 

As a first implication we summarise some consequences of the tightness estimates from ^3.51 more 
precisely of Corollary 13.241 and Lemma 13.251 

Remark 3.29. There exist constants d and D which depend only on the initial data g and T but not 
on no such that for all t G [0, T] we have 

M''°\t)>d, F("")(t)<D, M^'^'^°\t,a) < Dexp{-da). 

Our second result provides compactness in the space of bounded and continuous functions. 
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Lemma 3.30. For each n the set {fn ■ uq > is equi- continuous in BC{[Q, T] x [0, cxd)). 

Proof. Within this proof let e > be fixed and suppose that 6 is sufficiently small. We start with the 
modulus of continuity in a -direction as it can be estimated independently of n. Since the auxiliary 
Problem l3.8l is linear in / and invariant under shifts w.r.t a, we immediately obtain 

t 

S+f("o){t) - J^''»\t) = 5+5^"°) - + J r("«)(,s)r(t - ,s)j(5+/("")Gs) - ds, 



and Corollary 13.131 provides 

WSl&^Ht) - < ||5+.g("°) - < e. 

Therefore, all function fn^°\t) arc cqui-continuous in a-dircction (uniformly in n, hq and t). Now let n 
be fixed and consider uq > n. The first part of this proof implies 

\\{nS)f(-o)(t))^~fi-"Ht)\\^<e 

for all < (5 < (5. Note that here 5 is expected to depend also on n as the transport velocity increases 
with n. Moreover, by construction we have 

5 

+ =T(^)/("'')(i) + J r("'')(t + s)T((5-s)J/("°)(t + s)ds, 



and since r("f'(<) is uniformly bounded, see Remark 13.291 we infer that 

5 

||/("«)(t + <5)-T(5)/("°)(i)||^<C J \\{J.&°\t + s))J\^ds<C^S, 



with constant C„ depending on n. Finally, we have shown that 

\\fi""Ht + s)-&"Ht)\L<e 

for all sufficiently small S, and the proof is complete. □ 

We are now able to prove our main result, which was already stated in ^ Theorem l2.3l 
Lemma 3.31. There exists a mild solution f ^ y to Prohlem \2.1[ which has the following properties. 

1. f is non-negative with < \\g\\fj 

2. all states f{t) are admissible, so that P(f{t)) = 0, 

3. all states are regular in the sense of Definition \S.Sl\ 

4. f conserves the area with non-increasing number of grains. 

Proof. Let — [0, T] x [0, 00). The Arzela-Ascoli theorem implies that any equi-continuous subset of 
BC{Q) has a subsequence that converges to a bounded and continuous limit function on Q, and that this 
convergence is uniform on each compact subset of Q. Thanks to Lemma l3.30l there exists a sequence {nk)f. 
with Uk ^ GO such that f^''^ converges to a limit /2 locally uniform in BC{Q). Moreover, passing to a 
suitably chosen subsequence, still denoted by (nfc)j,, we can assume that 73"'''' fa for some /a S BC{^). 
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Iterating this argument, and using the usual diagonal trick, we finally find a subsequence along with limit 
functions / = (/2, ...,/„, ...) such that 

locally uniform in BC{n). 

for all n > 2. By construction, the limit / is non-negative and satisfies < HsH/j for all t. The 

uniform tightness estimates from Remark 13.291 imply the L^-regularity of /, and hence / G 3^, as well 
as M^{f{t)) < D exp {~dt). Moreover, since these tightness estimates control the tail behaviour with 
respect to n and a we infer that 

iv(/("'=)(<)) Nifit)), r(/("'«)(t)) r(/(t)), 

and 

p(/(n.)(t)) p(/(t)), A(/("'')(t)) Nifit)). 

From this we conclude that / is in fact a mild solution to Problem 12.11 a-nd the remaining assertions 
concerning A{f{t)), N{f{t)) and P{f{t)) follow from the corresponding properties of □ 

Remarks on well-posedness 

To establish well-posedness for our kinetic model we must prove both, the uniqueness of solutions and 
the continuous dependence on the initial data. This can done by means of energy methods, if the 'energy 
distance' of two states / and / is defined via 

r°° 1 1 

£ifJ) = Y.n exp(-a)(/„(a)-/„(a)) da + (7V(/) - iV(/)) . 

n>2 

The main result on £, which in turn implies uniqueness and continuous dependence for solutions to 
Problem [231 can be stated as follows. 

Lemma 3.32. For all < T < oo and any two mild solutions f and f from Lemma \3.31\ we have 

£if{t),fit))<expiCt)Eigr9) (50) 

for all < t < T , where the constant C depends only on T and the initial data g = /(O), g ~ f{0). 

The proof of ([50]) relies on careful estimates for (/(t), /(t)). Since a concise presentation would 
involve lengthy computations, we omit the details here and refer the reader to [Hen07j . 
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